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Abstract 

We provide a combinatorial description of exclusion statistics in terms of minimal 
difference p partitions. We compute the probability distribution of the number of 
parts in a random minimal p partition. It is shown that the bosonic point p = is 
a repulsive fixed point for which the limiting distribution has a Gumbel form. For 
all positive p the distribution is shown to be Gaussian. 
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1 Introduction 

Integer partition problem has a long history going back to Euler. The classical question 
asks: in how many ways p{E) can one partition an integer E into nonzero integer parts 
E = J2j such that hj > hj+i for all j = 1,2,.. .? For example, p(4) =5: 4 = 4 = 
3 + 1 = 2 + 2 = 2 + 1 + 1 = 1 + 1 + 1 + 1. Pictorially, one can represent the part hj as the 
height of the j-th column with nonincreasing heights such that the total height under the 
columns is E. Hardy and Ramanujan proved P that for large E, p{E) ~ 1 31/2^ 6°^, with 
a = 'K\f2j^. Similarly, one can ask the number of ways of partitioning the integer E into 
distinct integer summands, i.e., E = ^jhj such that hj > hj^i with strictly decreasing 
heigth. For example, the integer 4 can be partitioned into distinct summands in only two 
ways, 4 = 4 = 3 + 1. In this restricted case it is known that asymptotically for large E [2], 
p{E) ~ \-,j^e'^ where b = a/V2 = n/VS. 

Another way of representing integer partitions makes clear the connection with a gas 
of noninteracting quantum particles. Let be the number of columns of height h = i in 
a given partition, i.e. the number of times the summand i appears in a given partition. 
For example, in the partition 4 = 2 + 1 + 1, one has rii = 2, n2 = 1 and rij = for all 
j > 2. Then, E = J^i ^i^i where = i for i = 1, 2, . . . represent equidistant single particle 



energy levels and rij = 0, 1, 2, . . . represents the occupation number of the i-th. level. In 
the unrestricted problem, the occupation number rij = 0, 1, 2 . . . (Bosons) whereas in the 
restricted problem = 0, 1 (Fermions). Therefore, E = njej is the total energy of the 
system and 

p{E) = J2s[E-Y.n,eA. (1) 

rii \ 1=1 J 

If in addition, one restricts the number of summands to be A^, then the number p{E, N) 
of ways of partitioning E into parts is simply the micro-canonical partition function of 
a gas of quantum particles with total energy E and total number of particles N 

p{E,N) = Y,^(E-Y,n,tA WiV-f^nJ. (2) 

rii \ i=l J \ 1=1 / 

Evidently, p{E) = Yl^=QPiE,N). Even though the sum p{E) has similar asymptotic 
behavior for large E for Bosons and Fermions, i.e., \n{p{E)) ~ ^/E (up to a constant 
prefactor), we will show in this paper that p{E, N), as a. function of N for a fixed E, has 
rather different behavior for Bosons and Fermions. 

Thus a gas of non-interacting Bosons or Fermions occupying a single particle equidis- 
tant spectrum (ej = i) both have a combinatorial interpretation in terms of partitions of 
an integer E into parts. 

• Bose statistics corresponds to the case of unrestricted partitions rij = 0, 1, 2 . . .. 

• Fermi statistics corresponds to the case of restricted partitions with distinct sum- 
mands rii = 0,1. 

A natural question, that we address in this paper, is how to provide a combinatorial 
description of a quantum gas obeying exclusion statistics. Exclusion statistics is a gener- 
alization of Bose and Fermi statistics [H El El [7] . It has been found explicitly in quantum 
models of interacting particle systems, notably in the two dimensional lowest-Landau- 
level (LLL) anyon model [3] (i.e. the anyon model projected into the LLL of a strong 
magnetic field) and the one dimensional Calogero model [HI El [IHl IHl [121 US]- Note that 
the Calogero model can be obtained as a particular limit of the LLL-anyon model [1^ , the 
latter being a particular exactly solvable projection of the anyon model: it follows that 
exclusion statistics is deeply rooted in the more general concept of anyon statistics [15j . 
Unlike the Bose and Fermi statistics which describes noninteracting particles, a combi- 
natorial description of exclusion statistics is a priori quite nontrivial since the underlying 
physical models with exclusion statistics describe truly interacting A^-body systems. 

We show in this paper that a combinatorial interpretation of exclusion statistics in- 
volves a generalization of the partition problem known as the minimal difference partition 
(MDP) problem. In MDP, one partitions a positive integer E into A^ nonzero parts, 
E = J2f=i (with hj > for all j = 1, 2, . . . , A^) such that each summand exceeds the 
next by at least an integer p, i.e. [hj — hj+i) > p for all j = 1, 2, . . . , A^ — 1. Therefore, 
p = corresponds to unrestricted partitions (Bosons) and p = 1 to restricted partitions 
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(Fermions) into distinct parts. Even though the parameter p in MDP is an integer, one 
can analytically continue the results to noninteger values of p and we will show that 
for < p < 1, the MDP corresponds to a gas of quantum particles obeying exclusion 
statistics. 

Apart from establishing this equivalence between the MDP problem and exclusion 
statistics, we also provide a detailed analysis of the asymptotic behavior of Pp{E, N), i.e., 
the number of ways the integer E can be partitioned into N parts in the MDP problem, for 
all p > 0. This analysis tells us how the variable fluctuates from one partition to another 
for fixed E. Indeed, defining Pp{E) = J2n Ppi^^ total number of partitions of 

E and treating all such partitions equally likely, the ratio Pp{N\E) = pp{E , N) / pp{E) 
is the probability distribution of the random variable A^, given E. We show that this 
distribution, properly centered and scaled, has rather different limiting shapes for p = 
and p > 0. While for p = the scaled distribution is asymmetric and has a Gumbel 
shape, for p > (including the Fermionic case p = 1) the scaled distribution is symmetric 
and has a Gaussian shape. 

At this point, it may be useful to summarize our main mathematical results for the 
asymptotic behavior of Pp{N\E). For the Bosonic case {p = 0), the limiting shape of 
the distribution was first derived by Erdos and Lehner using rigorous methods involving 
upper and lower bounds [16j. In this paper we calculate the limiting shapes of Pp{N\E) 
for all p > 0. Moreover, our method allows us to compute the probabilities of atypical 
large fluctuations which go beyond the range of validity of the limiting distributions. 

For p = we show that Po{N\E), as a function of A^ for fixed E, has a peak at a 
characteristic value N^{E) ~ ^ a/E log(4£'/a^) for large E, where a = 71^/2/3, and the 
random variable A^ typically fluctuates around Nq(E) over a scale ~ ^/E. Moreover, in 
the vicinity of A'o(^) over a range |A^ - N^{E)\ ~ 0{VE), the distribution Po(A^I^) has 
a scaling form (or a limiting law). In terms of the cumulative probability, 

Qo{N\E) = J2 Po{N\E) ^fJ-^{N- N*{E)) , (3) 

where the scaling function Fq{z) has an asymmetric Gumbel form, thus recovering the 
Erdos-Lehner result [16j 

Fo{z) =exp[-exp[-z]]. (4) 

In contrast, for p > 0, the distribution Pp{N\E) has quite a different asymptotic 
behavior. It has a peak at a characteristic value N*{E) ~ ai{p)\fE and A^ typically 
fluctuates around N*(E) over a scale of ~ E^^^ for all p. Moreover we show that, on 
this scale, the fluctuations are Gaussian. More precisely, we show that in the vicinity of 
N*{E), the cumulative probability Qp{N\E) has a scaling form 

is a universal scaling function independent of p (> 0). The two nonuniversal scale factors 
ai{p) and a2{p) however depend explicitly on p and can be computed exactly. For example. 
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for Fermions {p = 1), we recover the Erdos-Lehner result for the mean 



ai(l) = ^ln(2) 



TT 



and get a new result for the variance 



a2(l) 



Stt^ - 36 ln^(2) 



.1/2 



(6) 



(7) 



Thus, as far as the limiting shape of the scaled distribution of Pp{N\E) is concerned, 
it is a universal Gaussian for all p > 0. The Fermionic case p = 1 is thus a representative 
of all p > and can be considered as an attractive fixed point along the p-axis (see Fig. 
([T]). In contrast, the Bosonic case p = represents a repulsive fixed point where the shape 
is Gumbel. 



p=0 p=l 

BOSON FERMION 

( GUMBEL ) ( GAUSSIAN ) 



Figure 1: Schematic flows along the p axis. The p = represents the Bosonic fixed point 
where the limiting distribution of Pp{N\E) is Gumbel. In contrast, the behavior for all 
p > is controlled by the Fermionic fixed point at p = 1 where the limiting distribution 
is Gaussian. 

The limit laws above describe the probabilities of typical fluctuations of around 
its characteristic value N*{E)7 In this paper, we have also investigated the probability 
of atypical large fluctuations of away from N*{E) and calculated the corresponding 
large deviation functions exactly. Like the limit laws, the large deviation properties for 
p > turns out to be rather different from the p = case, thus confirming the fixed point 
picture of Fig. [H Curiously though, we show that the large deviation function for any 
p > is related to that of p = via an exact nonlinear relation. 

The paper is organized as follows. In Section 2, we precisely define the MDP prob- 
lem, provide an exact derivation of the generating function of Pp{E, N) and establish a 
nonlinear relation between pp{E,N) with p > and po{E,N). In Section 3, we show 
how the MDP problem with < p < 1 corresponds to exclusion statistics. In Section 4, 
we provide detailed asymptotic analysis of Pp{E, N) for all p > and obtain the limiting 
shapes of the scaled distribution Pp{N\E) and also calculate exactly the associated large 
deviation functions. Finally, we conclude with a summary and open problems in Section 
5. 
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2 Minimal Difference Partition Problem: A Combi- 



natorial Approach 

In the minimal difference partition (MDP) problem, one partitions an integer E into 
nonzero parts, E = XljLi (with hj > for each j = 1,2,... A^) such that each part 
exceeds the next one by at least an integer p, i.e., {hj — /ij+i) > p for all j = 1, 2, . . . , — 1 
(see Fig. [2]). Let pp{E,N) denote the number of ways one can achieve this. Clearly, the 
cases p = and p = 1 reduce respectively to the unrestricted partitions (Bosons) and the 
restricted partitions (Fermions). The generating function for pp(E,N) is well known [3] 
and is given in Eq. f|TT]) . However, here we provide a simple derivation of this result that 
brings out in a direct way a nontrivial connection between the cases p > and p = 
which will be used later for the analysis of the asymptotic behavior. 



■ 5iP 



E= Ih. 

j=i ^ 



}>0 



Figure 2: A typical partition configuration of the MDP problem with N = 5. The column 
heights hj > for all j = 1, 2, ... A^ and their total height is E = Ylf=i ^j- addition, 
they satisfy the constraint, {hj — /ij+i) > p for an integer p for all j = 1, 2, . . . , A^ — 1. 

Let us first establish an exact one to one correspondence between a partition con- 
figuration of the MDP with nonzero p > and a partition configuration with p = 0. 
Let {hj} denote the set of nonzero heights in the partition of E = X^^i p = 

(Bosonic case). Thus, hj > hj^i for all j = 1, 2, . . . , A^ — 1. Let us now define a new set 
of heights h'j = hj + p{N — j) for j = 1,2, . . . , N. Thus, h'j — h'-j^^ = hj — hj^i + p for 
all J = 1, 2, . . . , A^ — 1 and /i'^ = /iat > 0. The new heights thus satisfy the constraint 
{h'j — h'j^i) > p for all j = 1, 2, . . . , A^ — 1 and their total height is given by 

N N 

E' = ^h'j = E + pN{N - l)/2 = ^hj+ pN{N - l)/2. (8) 
j=i i=i 

Therefore, the primed heights correspond to a partition configuration of the integer E' 
into A^ parts with p > 0. This exact correspondence then provides us with the following 
identity valid for all A^ 

Pp{E, N) = po{e-^- N{N -1),n). (9) 



5 



Thus if one can compute the partition function pq{E,N) for the Bosonic case {p = 0), 
this identity can be used to obtain exact results for any arbitrary p > 0, including the 
Fermionic case p = 1. 

For the Bosonic {p = 0) straightforward calculation gives the generating func- 

tion 

oo N 

Using the correspondence in Eq. Q gives the general result for all p > 

^N+pN{N-l)/2 



°° N+pN{N-l)/2 
E=l 



It turns out to be convenient sometimes to use the cumulative partition function 
Cp{E,N) = "^^i^Q PpiE, N). For p = 0, its generating function can be easily derived 
from Eq. (fTOj) and has a particularly simple form which turns out to be rather useful, 

CXD ^ 

g^°<^-^'-'- (l-.)(l-.^)...(l-:^") - '''' 
Comparing Eqs. (fTTj) and (|T^ one gets another identity 

Pp(E, iV) = Co (i? - iV - I iV(iV - 1), iv) (13) 

which we will be using later. 



3 The MDP with < p < 1 and Exclusion Statistics 

In this section we show that the MDP problem with integer parameter p, continued an- 
alytically to the range < p < 1 corresponds to a quantum gas of interacting particles 
obeying exclusion statistics. This correspondence is established at two levels: (i) at a 
microscopic level where we show in subsection 3.1 that the pp{E, N) of the MDP problem 
corresponds precisely to the micro-canonical partition function of the one dimensional 
Calogero model in an external harmonic potential and (ii) at a more general thermody- 
namical level in section 3.2. 



3.1 Equivalence between the MDP problem and the spectrum 
of the Calogero model in a harmonic well 

The aim of this subsection is to show that there is an exact one to one correspondence 
between the energy levels of the one dimensional Calogero model in an external harmonic 
well and the partition configurations of the MDP problem, continued analytically to 
< p < 1 in the sense explained below. 
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The Calogero model (for a review see [17|, [18J) describes an interacting quantum 
particle system on a line where the particles attract each other by an inverse square 
potential. In order to have a proper thermodynamic limit, one can either put the 
particles in a finite box of size L and then take the N ^ oo, L oo limit keeping the 
density N/L fixed. Alternatively, one can keep the particles on the infinite line, but switch 
on an external harmonic potential of strength u. In the latter case, one has to eventually 
take the limit — >■ in a suitable way. It turns out that while the model in a box is not 
integrable, the model in a harmonic potential is integrable. Setting the Planck's constant 
h = 1 and the mass of each particle m = 1, the quantum Hamiltonian of the model is 



JV^ ^2 . . TV 

H 



i= * i<j ^ * ^' 1=1 

where Xj represents the position of the i-th particle, a; represents the frequency of the 
external harmonic well and a G [—1,0] represents the coupling strength of mutually 
attractive interaction between the particles. In addition, the many body wavefunction 
must vanish at Xi = Xj for any pair {i 7^ j) of coordinates for a ^ 0. It turns out that 
while a = represents non-interacting Bosons (where the many body wavefunction is 
symmetric under the exchange of Xj and Xj), a = —1 represents noninteracting Fermions 
(where the wavefunction vanishes at Xi = Xj). For other values of a, this model is known 
to exhibit fractional statistics (see in particular in the next subsection its manifestation 
in the thermodynamics of the model). 

The many-body energy spectrum of this model is known exactly [8]. The energy 
E{{hj}) is labelled by non-increasing integers hi >h>h...>i 



^ 1 

J2h,--aN{N-l) 



2 



(15) 



By making a shift as in Eq. ([8]), i.e., defining a new set of variables h'j = hj + a{N — j), 

one can express the energy as E = u XljLi ^'j with the constraint that {h'j — /ij+i) > — «. 
Thus the spectrum of the Calogero model in a harmonic potential corresponds exactly to 
partition configurations of the MDP with parameter p = —a, but now p is a real number 
such that < p < 1. Hence, the micro-canonical partition function of the Calogero model 
pcai{E, N) denoting the number of configurations with energy E and number of particles 
N is directly related to the number of partitions Pp{E, N) of the MDP model via 

PcAE,N)=pp{E/uj,N) (16) 

This implies that the grand-canonical partition functions of the two models are also re- 
lated. Let Zcai(/9, z) = ^ pca.i{E, N)e~^^ be the grand-canonical partition function 
in the Calogero model in a harmonic well of frequency cj, where j3 is the inverse temper- 
ature and z is the fugacity. Similarly, we define Zp[(5, z) = J2e n Ppi^, N)e~^^ z'^ as the 
double generating function in the MDP problem with parameter p. The relation in Eq. 
(1161) then translates into the following relation between the grand partition functions 

Zc.i{l3,z) = Zp{ul3,z). (17) 
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3.2 Thermodynamic equivalence to exclusion statistics 

Exclusion statistics can be most conveniently defined in the following thermodynamical 
sense. Let Z{P, z) denote the grand partition function of a quantum gas of particles at 
inverse temperature (3 and fugacity z. Such a gas is said to obey exclusion statistics with 
parameter < p < 1 if z) can be expressed as an integral representation 



where p(e) denotes an effective single particle density of states and the function yp{x), 
which encodes fractional statistics, is given by the solution of the functional equation 



Note that for p = 0, one gets yp{x) = 1/(1 — x) and for p = 1, yp{x) = (1 + x). In 
these two extreme cases, Eq. (fTSll reduces to the standard grand partition functions of 
noninteracting Bosons and Fermions respectively. The fractional statistics with parameter 
< j9 < 1 (that corresponds to an interacting gas) then smoothly interpolates between 
these two extreme cases. 

There are at least two microscopic quantum models whose grand-canonical functions 
have the form of Eq. f|T8l) . The first example is the LLL anyon model [5| in the infinite 
volume limit which can be shown to satisfy Eq. ffTHj) with an effective density of states 
p(e) = 6{e — uJc) where B is the external magnetic field, 0o = 27r/e is the flux quantum, 
ujc = eB/2m is the cyclotron frequency and V is the infinite area of the system. In this 
model, the parameter p = corresponds to the the flux carried by each anyon in units 
of the flux quantum. The second example corresponds to the one dimensional Calogero 
model defined in Eq. (fill) again in the infinite box limit. In this case, one can show that the 
grand partition function again can be written in the form as in Eq. ffT^ with an effective 
single particle density p(e) = L/VSn'^e where L is the infinite length of the system. In 
both cases, the thermodynamics is computed in the presence of a long distance harmonic 
well regulator, and the thermodynamic limit where the external frequency ^ is taken 
in such a way so that one correctly recovers the infinite box limit. 

Here we show, using the equivalence to the MDP problem in Eq. 0171) . that the grand 
partition function of the one dimensional Calogero model in an external harmonic well of 
frequency u, in the limit u; — *■ can again be written in the general form as in Eq. (ITSll . 
but now with an effective constant density of states p(e) = l/u. Note that this is different 
from the Calogero model in a infinite box of size L (the second example mentioned in the 
preceding paragraph): here, the particles are sitting inside a harmonic well with almost 
vanishing but non zero frequency. 

To proceed, we first calculate the grand partition function of the MDP problem, 
Zp{iu(3,z) = T.E,N PiAE,N)e-^^'^z^ , starting from Eq. We set x = e'^^ in Eq. 

(flTi) . multiply it by z^ and sum over N . Next we take the logarithm on both sides and 
then make a cluster expansion, \n Zp{uP, z) = Yl'^=i^nZ"' ■ Now, taking the u; — > limit 




(18) 



yp{x) -xyl P(x) = 1. 



(19) 
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(keeping (3 fixed), one gets 

&.>2-^^n(i-f) (20) 

up up k 

Consequently, since from (fT^ \nyp{x) = J2'^=i IT 11^=1(1 ~ one infers that, provided 
the series is convergent, that is ze''^^ < 1 

/oo 
In yp{ze-'^'^')de. (21) 

Making a further change of variable ue —>■ e, it follows that in the limit u 

1 f°° 

In Zp{u(3, z) ^ - lnyp{ze-'^')de. (22) 
^ Jo 

This is again of the form in Eq. (|T8l) with p(e) = 1/u. Using the equivalence in Eq. ( ITTj) 
we then conclude that the Calogero model in an external harmonic well with vanishing 
frequency, which precisely corresponds to the MDP problem with parameter < p < 1, 
can be viewed as a gas of particles obeying exclusion statistics with a statistical parameter 
a = —p and a constant density of states. 



4 Partition Asymptotics in MDP with p > 

In this section we compute explicitly the asymptotics of the probability distribution 
Pp{N\E) in the MDP problem for all p > 0. We show that while the limiting shape 
of this distribution (properly centered and scaled) is Gumbel for p = 0, it is Gaussian for 
all p > including the Fermi case p = 1. 

4.1 Bosonic case p = 

Our starting point is the generating function for the Co{E,N) in Eq. ( fT2l) . We formally 
invert this generating function using Cauchy's theorem and write 



2m J x^+^ (1 - x)(l - a;2) . . . (1 - x^) 

N 

/3E-^ln(l-e-'^'=) 



2^ / 



k=l 



(23) 



where the integration is in the complex x plane along a contour around the origin and we 
have made a change of variable x = exp(— /?) in going to the second line. For large E, one 
can then analyze the leading asymptotic behavior by employing the saddle point method 
in the complex /3 plane. Anticipating that for large E, the most important contribution 
to the integral will come from small (3, we first obtain the leading small /3 behavior of 
the action Se,n{,P) = PE — J2k=i^^ ~ e"'^'^). Using the Euler-Mclaurin summation 
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formula, one can easily show that in the limit of (3 0, N ^ oo limit but keeping jSN 
fixed, the action can be written as 



1 r'^^ tdt 

SeMP) c:, f3E + - / -iVln(l 
P Jo e — i 



-0N 



We next maximize the action with respect to (3, i.e., we set dS/dP = to get 

tdt 



For a given large E, one gets P* by implicitly solving the saddle point equation 
substitute it back in the action Se,n{P*)- Thus, to leading order, 

Co(^,iV) ~exp [SeAP*)] 
where the micro-canonical entropy Se,n{P*) can be written as 



13* N 



tdt 
e* - 1 



(24) 

(25) 
and 

(26) 
(27) 



To bring out the scaling form of Cq{E, N) explicitly for large E and A^, we next proceed 
as follows. It is evident from the structure of the saddle point solution that P* ~ 
for large E, whereas P* ~ for large indicating that the correct scaling variable is 
X = N/ \/E. Next we set P*N = H{x). In terms of these new scaling variables, the saddle 
point solution in Eq. fl25l) can be recast as 



H^ix] 



x^ 



tdt 



{2i 



Thus, given x, one has to find H{x) by implicitly solving Eq. ( l28l) . The entropy in Eq. 
(1271) becomes Se,n{P*) = ^/E g{x) where the scaling function g{x) is given from Eq. fl27j) 

as 

^(a;) = 2:^-a;ln(l-e-^(")). (29) 

Thus, asymptotically for large and E, keeping the ratio x = N/ \fE fixed, the cumula- 
tive number of configurations Co{E, N) for Bosons can be written as 



Co(E,Ar) ~exp 



Eg 



N 



(30) 



where g{x) is the large deviation function given exactly by Eqs. (l29l) and (l28l) . This is 
the main result of this subsection. 

The function g{x) has to be determined numerically by solving the implicit equations 
fl^^ and fpH|) . A plot of this function is given in Fig. The asymptotic properties of 
g{x) for small and large x can be worked out easily. It can be shown that 



g{x) ~ — 2xln(x) as a; ^ 
2 

a exp(— ax/2) as x 



oo 



(31) 
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Figure 3: The large deviation function g{x) for Bosons {p = 0). 



where a = ny^2/3 = 2.5651 .... 

The result g{x oo) = a implies, from Eq. ( |30l) that p{E) = C{E,N — > oo) ~ 
exp[a\/E] to leading order for large E, thus recovering the famous Hardy- Ramanuj an 
result [T]. The normalized cumulative distribution of (given E), Qq{N\E) = Co{E, N)/p{E) 
then has the large deviation form 



Qo(iV|^)-exp 



-VE(I> 



N 



where = a — g{x) 



(32) 



and ^(x) has the asymptotic behavior 

$(x) ^ a + 2xln(x) 



- exp{—ax/2) as x 



oo. 



(33) 



As a; oo, i.e., as >> clearly Q{N\E) ^ 1 as expected, since it is the 

normalized cumulative distribution of A^. The precise approach to 1 can be obtained 
using the large x asymptotics of $(x) in Eq. ( !33l) . Substituting this behavior in Eq. ( l32l) 
one gets for >> ^/E, 



QoiN\E) ~ exp 



E exp{-aN/2VE) 



2y/E 



{N-N*iE)) 



(34) 



where the characteristic value of the random variable A^ is Nq{E) ~ ^yE log{4:E/a'^) 
and the scaling function has the Gumbel form, -Fo(^) = exp[— exp[— z]]. Evidently, the 
probability distribution PoiN\E) = QiN\E) - Q{N - 1\E) ^ dQoiN\E)/dN has the 
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scaling form 

PoiN\E) ^ (N - N*iE))^ where = exp[-;. - exp[-z]] (35) 

which is highly asymmetric around the peak at = Nq{E). This limiting distribution 
of that describes the probability of typical fluctuations of of ~ 0{\fE) around the 
peak at Nq{E), was originally derived by Erdos and Lehner by computing upper and lower 
bounds to the probability [15]. Our method allows us to obtain a more general result in 
Eq. (1321) which is valid over a wider range and reduces to the Gumbel limiting form near 
the peak. A rigorous mathematical derivation of this result, including the exponential 
prefactor, can be found in the work of Szekeres [T9] . 



4.2 The case p> 



For p > 0, one can directly obtain the asymptotic behavior of pp{E,N) by using the 
identity in Eq. ( |T3|) and the already derived asymptotic behavior of Co{E,N) in Eq. 
(l30l) . In the scaling limit when A^ and E are both large but the ratio x = N/ ^J~E is kept 
fixed, one gets to leading order 



Pp(E, A^) ~ exp 



N 



with 



fp{x) = a/1 -px'^/2 g 



X 



y^l - pX^/2 



(36) 



(37) 



This is 



where g{x) is given in Eqs. ( l29l) and ( l28l) . 

Note that the function fp{x) has nonzero support only over x G 0, \/2jp 

easy to understand from the fact that for p > 0, has a minimum value for any given 
A^, or equivalently A^ has a finite maximum value for any given E. For example, for the 
Fermionic case (p = 1), the lowest value of E for a given A^ corresponds to the Fermi 
energy Ep = N{N + l)/2 where one puts one Fermion at each single particle level ei = i 
for 2 = 1,2, ... ,A^. Thus, E > N{N + l)/2 for all A^. In other words, for large A^, 
A^ < V2E, i.e., X < \f2. Similar arguments can be given for any positive p > 0. Unlike 
the function g[x) which is monotonically increasing, the function /p(x) in Eq. (1371) is a 

non-monotonic function in a; G 



It vanishes at the two ends as 



/p(x) ~ —1x ln(x) as x 







TT 



V6, 



(2p)vy^- 



X 



(38) 



and has a unique maximum at x*(p) = ai(p), where ai(p) can be obtained by setting 
dfp{x)/dx = in Eq. (1H7|) and then using the known properties of g{x). A plot of fp{x) 
for p = 1 (Fermi case) is given in Fig. (jll). 
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Figure 4: The large deviation function fi{x) for Fermions {p = 1). 



By playing around with the form of fp{x) in Eq. fl37p and that of g{x) in Eqs. 
and (l28l) one can derive a number of explicit results. We skip the details here and just 
mention the results. For example, the location of the maximum x*{p) = ai{p) is given by 

ai{p)= ^""^^ (39) 

v/p(ln2/o^)V2 + ^^2(l-l/yS) 

where — y^~^ = 1 and Li2{z) = Ylh=o dilogarithm function. For example, 

in the Fermionic case p = 1, we get = 2 and ai(l) = 2v/31n(2)/7r = 0.764304.... 
Similarly, the value of the function at the maximum fp{x = ai{p)) can be shown to be 



fp{x = ai(p)) = 2 v/i:^2(l-l/z/S)+p(ln 1/0^)2/2 (40) 

For example, for p = 1, it gives fp{x = ai(l)) = vr/v^. For an arbitrary p, this formula 
which goes back to Meinardus [2T1| , provides a generalization of the Hardy- Ramanuj an 
formula for p{E). One can check that Eq. fHCTl ) coincides with the integral formula obtained 
by Blencowe et al [20! '^^^ ^^^^ investigated the Haldane statistics but didn't provide any 
combinatorial interpretation of their result. Around the maximum at x = ai(p), the 
function fp{x) can be expanded in a Taylor series and up to the quadratic order 

fp{x) ~ /p(ai(p)) + (x - ai(p))2 + . . . (41) 

where 02 (p) can also be evaluated. For example, for p = 1, we get 



«2(1) 



37r2 - 36 W{2) 



0.478815... (42) 
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Evidently, one can easily evaluate the asymptotic behavior of pp{E) = "^j^f Pp{E, N) 
for large E by replacing the sum by an integral, use the large deviation form in Eq. flHUl) 
for Pp{E, N) and then using the saddle point method. To leading order, this gives 

Pp(E)~exp[/p(ai(p))v^]. (43) 

The normalized probability distribution of (for fixed E) Pp{N\E) = pp{E , N) / pp{E) 
then has the large deviation asymptotics 



Pp{N\E) ~ exp 



E^p 



N 



where ipp{x) = fp{ai{p)) - fp{x) (44) 



Thus for all p > 0, the probability distribution Pp{N\E) has a peak at a characteristic 
value N*{E) = ai(p)^/E (note the difference from the Boson case p = where Nq{E) ~ 
\/EId.(E)). Using the expansion in Eq. (HTj) it follows that in the vicinity of N*(E) (over 
a scale of ~ 0{E^^^)), Pp(N\E) has a Gaussian limiting form 



Pp{N\E) ^ — where F(z) = ^e'^^^ (45) 
Note in particular, that the standard deviation measuring the root mean square fluctua- 



tion of A^, (Jp{E) = y ((A^ — N*{E)y) grows with as a power law, ap{E) ~ a2{p)E^/^ 

where the exponent 1/4 is universal for all p > 0. Moreover, apart from nonuniversal p 
dependent scale factors such as ai{p) and a2(p), the full distribution Pp{N\E) also has 
the same universal Gaussian limiting form for all p > 0. Thus, the Fermi point p = 1 is a 
generic point that is representative of all values of p > as far as the limiting distribution 
is concerned. In this sense, all p > behavior is controlled by the attractive Fermi fixed 
point as shown in Fig. ([T]). The Bosonic fixed point at p = 0, on the other hand, is a 
repulsive one. 



5 Summary and open problems 

To summarize, in this paper we have provided a combinatorial interpretation of exclusion 
statistics in terms of minimal difference partitions(MDP). This correspondance is based 
on the observation that the grand-canonical partition function of the Calogero model 
coincides with the generating function of MDP. By going to the grand-canonical ensemble 
and taking a suitable thermodynamic limit, we have recovered the functional equation 
characteristic of exclusion statistics. Apart from establishing this correspondance, we 
have also provided a detailed analysis of the asymptotic behaviour of pp{E,N). Our 
approach uses a mapping with the bosonic problem which holds for arbitrary p G [0, 1]. 
In physical terms this generalises the well known mapping between fermions and bosons 
with a linear dispersion law [22l |23] . The fact that this mapping has a number theoretical 
interpretation was apparently not known before. By using this mapping, we obtain a 
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general description of the limiting laws of Pp{N\E) for all p > 0. We find that the bosonic 
point is a repulsive fixed point where the statistics is Gumbel. In contrast for all p > 0, 
the distribution is Gaussian. Several questions emerge from this work and would be worth 
investigating. 

1) The regime p < 0. In this case the functional equation ( fT9|) still holds where yp{x) 
can be interpreted as the generating function of connected clusters on a p-ary tree. A 
preliminary investigation of this model shows that the scaling behaviour of Pp{E, N) is 
quite different from the previous case [21]. 

2) In this work we have limited ourselves to the integer partition problem or equiva- 
lently to a quantum gas of particles with equidistant single particle spectrum, i.e., with 
a constant density of states p(e) = const.. It would be interesting to investigate general 
partitions of the form E = 'Y^Uii'^ that corresponds to having a power-law density of 
states, p(e) ~ e^^^~^. In this case we have shown in a recent work [2S] that the bosonic 
sector gives rise to the three universal distribution laws of extreme statistics, namely the 
Gumbel, Weibull and Frechet distributions. It would be interesting to explore the general 
p > case including the fermionic sector and see if the bosonic point p = is still a 
repulsive fixed point. 

3) For the bosonic case (p = 0), Vershik [2S] and Temperley [27] calculated the limiting 
shapes of the Young diagram, i.e., the average height profile for a fixed but large E. This 
result can be generalized [2lj to the case p > using the functional equation ( fT9l) . 



References 

G.H. Hardy and S. Ramanujan, Proc. London. Math. Soc. 17, 75 (1918). 



[2] 

[3] 

[4] 
[5] 

[6] 
[7] 



[9 



G.E. Andrews, The Theory of Partitions (Cambridge University Press, Cambridge, 
1998). 

G.E. Andrews, R. Askey, R. Roy, Special functions, Encyclopedia of Mathematics 
and its applications 71, (Cambridge University Press, Cambridge, 1999). 

F.D.M. Haldane, Phys. Rev. Lett. 67, 937 (1991). 

A. Dasnieres de Veigy, S. Ouvry, Phys. Rev. Lett. 72, 600 (1994); Modern Phys. 
Letters B 9, 271 (1995). 

Y.S. Wu, Phys. Rev. Lett. 73, 922 (1994). 

M.V.N. Murthy and R. Shankar, Phys. Rev. Lett. 73, 3331 (1994); ihid 72, 3629 
(1994). 

F. Calogero, J. Math. Phys. 10, 2191 (1969); ihid 10, 2197 (1969); ihid 12, 418 
(1971). 

J. Moser, Adv. Math. 16, 197 (1975). 



15 



[10] B. Sutherland, J. Math. Phys. 12, 246 (1971); ibtd 12, 251 (1971); Phys. Rev. A 4, 
2019 (1971). 



[11 
[12 

[13 

[14 
[15 
[16 
[17 
[18 
[19 
[20 
[21 
[22 
[23 
[24; 
[25 
[26 
[27 



S.B. Isakov, Int. J. mod. Phys. A9, 2563 (1994); ibid Mod. Phys. Lett. B8, 319 (1994). 
A.P. Polychronakos, Nud. Phys. B324, 597 (1989); Phys. Rev. Lett. 69, 703 (1992). 

F. D.M. Haldane, Z.N.C. Ha, J.C Talstra, D. Bernard, and V. Pasquier, Phys. Rev. 
Lett. 69, 2021 (1992). 

S. Ouvry, Phys. Lett. B510, 335 (2001) 

J. M. Leinaas and J. Myrheim, Nuovo Cimento B37, 1 (1977) 

P. Erdos and J. Lehner, Duke Math. J. 8, 335 (1951). 

V. Pasquier, jhep-th/9405104 (unpubhshed) . 

A.P. Polychronakos, J. Phys. A. 39, 12793 (2006). 

G. Szekeres, Quart. J. of Math. (Oxford)(2) 2, 85 (1951); tbid 4, 96 (1953). 
M.P. Blencowe and M.P. Koshnick, J. Math. Phys. 42 5713 (2001). 

G. Meinardus, Math. Zeitschr. Bd 61, 289 (1954). 
D.V. Anghel, J. Phys. A Math.Gen. 35, 7255 (2002). 
K.R. Patton, M.R. Geller, M.P. Blencowe 357, 427 (2005). 
A. Comtet, S.N. Majumdar and S. Ouvry, unpublished. 

A. Comtet, P. Leboeuf and S.N. Majumdar, Phys. Rev. Lett. 98, 88301 (2007). 
A.M. Vershik, Russ. Math. Surveys 52, 379 (1997). 

H. N.Y. Temperley, Proc. Roy .Soc. A 48, 683 (1952). 



16 



